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Abstract. - We consider the increase of the spatial variance of some inhomogeneous, non- 
equilibrium density (particles, energy, etc.) in a periodic quantum system of condensed matter- 
type. This is done for a certain class of initial quantum states which is supported by static 
linear response and typicality arguments. We directly relate the broadening to some current auto- 
correlation function at finite times. Our result is not limited to diffusive behavior, however, in 
that case it yields a generalized Einstein relation. These findings facilitate the approximation of 
diffusion constants/conductivities on the basis of current auto-correlation functions at finite times 
for finite systems. Pursuing this, we quantitatively confirm the magnetization diffusion constant 
in a spin chain which was recently found from non-equilibrium bath scenarios. 
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, Any current of some physical quantity like, e.g., elec- 
trons through solids is either induced by an external me- 
chanical force F (e.g., electric field) or a spatially non- 
uniform density (density gradient Vp). For systems fea- 
turing normal transport the currents in those two cases 
are routinely assumed to be determined by 



kF , 



(1) 



where j denotes the current; k the conductivity; and V 
the diffusion constant. As well-known, the derivation of 
the l.h.s. of Eq. (1) from linear response theory is rather 
Straightforward and leads to the Kubo formula [1] which 
nowadays is a standard approach to transport in quantum 
systems [2-7]. However, a direct derivation of the r.h.s. of 
Eq. (1) from the corresponding scenario, i.e., without any 
external force appears to be more subtle. The overview 
paper [8] by Zwanzig, e.g., lists essentially six families of 
approaches, each of them based on different sorts of as- 
sumptions such as Fokker-Planck dynamics in the space of 
relevant observables [9-11], Onsager's regression hypothe- 
sis [12,13], restriction of the dynamics to local equilibrium 
states [14], etc. These subtleties are especially disturbing 
in the context of heat conduction, since there simply is 
no practical external force which could cause a heat cur- 
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rent [15-19]. However, all of the above approaches even- 
tually give the diffusion constant in terms of current auto- 
correlation functions also, i.e., proportional to the conduc- 
tivity. Of course, this is what Einstein and Smolouchowski 
firstly suggested in their ground-breaking work on Brow- 
nian motion. The above approaches are comprehensively 
covered in textbooks like, e.g., [1]. Therein the interested 
reader may also find a discussion of their consistency and 
implications. In order to motivate our present alterna- 
tive approach to this extensively debated subject, we sim- 
ply discuss here a standard formulation rather than going 
through the above discussion. In [1], e.g., one finds the 
following expression for a particle diffusion coefhcient 



kTdn/d^ 



dte"^' lim 

g-.0 
L — 'OO 



(2) 



where n is the equilibrium particle density; ^ the chem- 
ical potential; J' denotes the "random current" ; its 
Fourier component with wavevector q; the brackets en- 
code the Kubo-inner product; and is the volume of 
the system. Conceptually, the most challenging part 
is probably to concisely show that in the above limit 
the auto-correlation of the random current may be re- 
placed by the auto-correlation of the "true current", i.e., 
{J'-qiQ); Jq{t)) {J-q{0); Jq{t)) (slucc the former in- 
volves the "irrelevant part" in terms of projection meth- 
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ods), however, this replacement appears to be generally 
accepted. Practically, the limits g — > 0, a; ^ imply that 
only the infinitely slow dynamics of density structures of 
infinite length scale can be expected to behave diffusively 
with the above diffusion coefficient. So. even if one knew 
the exact current auto-correlation function for the infinite 
system, that would not bare implications on the dynam- 
ics at finite time- or length-scales. And, even worth, if 
one has some information on the current auto-correlation 
function but only up to a finite time, (2) does not allow 
for any conclusions on the dynamics. Both issues become 
manifest, if the integral in (2) does not converge, then no 
information results except for the "non-diffusiveness" of 
the dynamics. Furthermore, it may be a little subtle to 
generalize (2) to, e.g., a micro canonical ensemble. 

Thus, in this paper we approach the subject neither 
from projection techniques nor from linear response. In- 
stead we (somewhat arbitrarily) "coarse-grain" the quan- 
tum system spatially into subunits. The density profile 
and the current operator are then discretely formulated 
on the basis of this coarse-grained description. After- 
wards a specific class of initial states featuring such a 
non-uniform density profile is introduced. Then, simply 
by applying Heisenberg's equation, the evolution of the 
variance of the density profile is expressed in terms of a 
double temporal integral of the current auto-correlation 
function. Based on this result, we discuss the connection 
between diffusion coefficient and conductivity. We further 
discuss the implications of our specific choice for the initial 
state w.r.t. external perturbations and quantum typical- 
ity. Finally, we numerically calculate the integral of the 
magnetization current auto-correlation function in a XXZ 
spin chain and quantitatively compare the outcome to re- 
cent results on the diffusion constants for such a system 
from non-equilibrium bath scenarios. 

The above mentioned periodic spatial coarse-graining 
scheme is introduced to facilitate a consistent definition 
of a local current. It is most conveniently explained for 
(but not limited to) an one-dimensional system. To those 
ends the transported quantity X as well as the Hamilto- 
nian H are decomposed into formally identical addends 
which correspond to different positions, i.e., ^^x^ = X, 

= H. Thus, is a local density of the trans- 
ported quantity. Note that the may be defined on 
or in between the positions of the x^, or both. We con- 
sider quantities which are conserved on the full system, 
i.e., [if,^] = 0. We further require 



= i[H, Xfi] = , Xfi] + t[h 
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Here, /i^- [h^+) is supposed to denote the local sub- 
unit of the Hamiltonian which is located directly on the 
l.h.s. (r.h.s.) of x^. This implies a kind of locality. How- 
ever, such a description may always be at least approxi- 
mately enforced, if the interactions are reasonably short- 



ranged. It may require the usage of ft.^ that are larger than 
a single elementary cell. Routinely, the comparison with a 
continuity equation suggests a definition of local currents 
according to the scheme indicated in Eq. (3) [2,3,6,20]. 
This is consistent, if , x^] -\- [h(^^j^iy , i (p+i)] = 0. The 
latter holds, if X is globally conserved. 

We now define the class of initial states we are going to 
consider. Those read 



(ip). (4) 



Here, Peq is any stationary (equilibrium) state of the full 
system, i.e., [i?,pcq] = 0. The brackets (...) denote full 
equilibrium averages, i.e., {A) = Tr{ylpoq}- Let {{A;B)) 

indicate the inner product {{A; B)) = Tr{A pcq B pcq}. We 
then define 

c{t, iJ-i^) = \{{d^it);d,)) , £2 = ^((4 (t); J,)) (5) 



which clarifies the from (4). does not depend on time, 
for periodic systems it further does not depend on ly. For 
interpretational reasons we note that in this case may be 
rewritten as £2 = {{X ^JXJ; X ~ {X))) / L = {{X-{X) f)/L 
(here we exploited [A",^f] = 0), where now and in the 
following L indicates the number of subunits in the full 
system, hence quantifies the equilibrium fluctuations 
of the transported quantity. Thus, w.r.t. (2), this implies 
e2 = kTdn/d^. 

We denote the actual expectation value of the local de- 
viation of the transported quantity from full equilibrium 
by d^{t), i.e., 

d,{t)^Tr{d,{t)pm- (6) 
This way we may write 

U fl l^L 

We are going to analyze the spatial variance of those de- 
viations from equilibrium, while we require them to be 
normalized to one, i.e., X^/j'^mC^) ~ ^- "^^^^ normalization 
is implemented by a corresponding choice of the Then 
we may directly, quantitatively compare to a discrete dif- 
fusion equation, as outlined below. The above mentioned 
spatial variance W'^{t) simply reads 



(8) 



If now, hypothetically, the dynamics of the (t) were gen- 
erated by a discrete diffusion equation of the form 

4(t) = V{t) [d^-i(t) -2d^{t) + d,,+i{t)], (9) 

then the evolution of the variance would read 

[mm^2vit) (10) 
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which holds, if the d^{t) vanish at the ends of a chain or 
are reasonably concentrated at a sector of a ring. (Such 
a concentration will be assumed throughout this paper.) 
Exploiting this, we arc able to deduce a diffusion constant 
X'(t) from the evolution of the variance. To those ends we 
rewrite the variance using (7) which yields 

- [J2^^5,]' ^[Y,vcit,fi)]\ (11) 

If the system features "space inversion symmetry" , which 
we assume in the following, the last term on the r.h.s. of 
the above Eq. (11) vanishes. In order to relate the evo- 
lution of W'^{t) to a current auto-correlation function, it 
turns out to be helpful to consider its second derivative 
w.r.t time. According to (11), this reads 

[Wm-T.'^'cit,rj). (12) 

We may evaluate this using Heiscnberg's equation: 

c{t,rj)^^^{{[H,[H,d,it)]];do)) (13) 
which may be rewritten as 

cit,ri) = -\ UH,d^{t)];t[Hjo])) (14) 
or, exploiting (3) and (4), as 

c(t, V) = iiOrj-i - 'in); 0-1 - 3o))) ■ (15) 

Inserting the above Eq. (15) into (12), and exploiting again 
that {(jri+((t)] jri)) does not depend on rj and further van- 
ishes for C i, we obtain 

[wm^^{{Jit);J))- (16) 

Here, J denotes not a local but the total current in the full 
system, i.e., J = X^^J/^- According to (10), the diffusion 
constant corresponds to the first derivative of the variance 
which reads 

[mm^[W^]it^0) + 2 f At'^{{J{t')-J)). (17) 

We suggest here to assume that [M^^J^^ = 0) = 0. This 
surely holds for T — > co, since in that limit W^{t) is sym- 
metric w.r.t. time. However, irrespective of T, according 
to typicality arguments, there are overwhelmingly more 
(pure) states corresponding to higher W'^ compared to 
any lower W"^ [21]. Thus, it is not to be expected that 
any state evolves towards lower , forwards or back- 
wards in time, unless it has been deliberately constructed 



to do so. We hence expect W'^[t) to be essentially sym- 
metric w.r.t. time, irrespective of T. Comparing to (10), 
we eventually conclude for the diffusion constant 

V{t)^ f At'^{{J{t')-J)), (18) 

which is the first main result of this work. 

In order to compare this to force-driven transport (as 
calculated from the Kubo formula) , we write out the inte- 
grand in (18) explicitly finding 

-^((J(t);J)) = V'(pmP„)5||J„„||2cos(u;,„„t) 

(19) 

with pm = (peq)mm, Jmn as thc matrix elements of the 
respective operators in the energy eigenbasis and with 
ojmn = ~ as the difference of the energy eigen- 
values En- Force-driven transport (within the linear 
regime) is routinely described by a response function that 
relates the current density to the external force: 

f dt'$(t')F(t'). (20) 

J —oo 

For a canonical equilibrium the low- frequency part of $(i), 
for which w„m ^ kT applies, may be written as 

$'°"'(t) ^ J^Y.^P^P^)'^\\J^n\? CO&{u0^nt) . (21) 

Thus, defining the response to slowly varying fields as 
a{t) = /g dt' $(i'), we get for times t > h/kT 

kT 

V{i)^—a{t) (22) 

which is a somewhat generalized Einstein relation, since 
it applies to all quantities (not only particles) . Obviously, 
in the case of diffusive transport the time dependencies of 
X'(t) and a{t) arc expected to vanish. 

In thc following we will discuss thc choice of thc initial 
non-equilibrium state in Eq. (4). This choice is essentially 
supported by two arguments. 

(I.) Assume the system was exposed (before any 
transport dynamics starts) to an, additionally weak, 
static "potential" such that the perturbed Hamiltonian 
reads H' = H + "^^v^x^. Then the corresponding 
canonical equilibrium state surely is given by p'^^ — 
cxp{~H'/kT)/TT{cxp{~H'/kT)}. (Note that this is of 
thc same form as a so-called "local equilibrium state"). 
Now, yOgq may be calculated for small using static lin- 
ear response [1]. Doing so, one finds that "low- frequency" 
contributions of thc "non-homogeneous" parts of p'^^ and 
p{0) are proportional to each other, i.e., 

{P'cq - Poq)m« OC (p(0) - Pcq)mn 

for \ujmn\ < kT , OC Sf^ (23) 
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with poq as the canonical state w.r.t. to H alone. All the 
dynamics of and p(0) under H stem from their inho- 
mogencous parts. All dynamics for times t ^ h/kT is 
controlled by their low-frequency contributions. In other 
words: if the initial density profile is induced by a previ- 
ous static external potential (which is then removed), the 
broadening of this density profile can be expected to be of 
the same form as the broadening corresponding to the hy- 
pothetical initial state p(0), as discussed in detail above. 
Our considerations therefore apply to this frequently dis- 
cussed type of initial state. 

(II.) A possibly even stronger motivation for the con- 
sideration of the initial states in Eq. (4) may be given 
along the lines of typicality [21-23]. Recently, it has been 
pointed out that in quantum systems with large Hilbert 
space dimensions the evolutions of observables of the type 
{^{t)\A\^{t)) do not depend crucially on the details of the 
initial state |V'(0)) [24]. Concretely, this means that the 
overwhelming majority of a set of initial states |'0„(O)) 
which all feature the same initial expectation value a of 
an observable A, i.e. (V'n(0)|A|'!/'n(0)} = a, will approxi- 
mately yield the same evolution of this expectation value, 
namely (V'„(i)|i|V'n(0) ~ Tr{ip™(^)}^ ^^l^^^ ^hE(^) 
represents a Hilbert space ensemble and explicitly reads 
p^^(O) = 1-1- JA, where / = /(a) is a pertinent scalar 
function of a. Although not proven so far, it is natural 
to suggest a generalization of this dynamical typicality 
to many observables such that for an initial set specified 
by (V'n(O)|Am|'0n(O)) ~ cLyn there exist typical evolutions 
for the expectation values of the type {^n{t)\Am\^n{t)) ~ 
Tr{i™pHE(^)| but now with a pHE(o) = 1 + ^„ i,„, 
where the /,„ — fm{{o-i}) are functions of the a;. If we 
consider for the moment a micro canonical ensemble in 
(4), we have to choose = Il£;/Tr{Il£;} with be- 
ing a projector projecting out the subspace spanned by 
the energy eigenstates corresponding to energies within 
a certain interval around E. With this choice the initial 
state in (4) reads p{0) = flsil + E^ S^,/e^ df,)IlE/Tr{IlE} 
which essentially is the same state as the latter p^^(O) 
with Afi = within the respective subspace. Thus, even 
if the true physical initial state may not be of the pre- 
cise form in (4), the initial state in (4) nevertheless gen- 
erates the typical evolution of the spatial density profile. 
Or to rephrase, if some initial state features a given spa- 
tial density profile and lives within a given energy region, 
the dynamical broadening of the density profile will most 
likely be described by (17), irrespective of the details of 
the initial state. 

Finally, we compare approximations to the diffusion 
constant 2?(i), as given by Eq. (18), with recent results 
in the literature. We approximate 2?(t) for finite times, 
simply by numerically exact diagonalization of finite sys- 
tems. Concretely, we analyze magnetization transport, 
i.e. = in the anisotropic s = 1/2 Heisenberg 




Fig. 1: The diffusion constant T>{t), as defined by Eq. (18), 
for magnetization transport (spin transport) in the anisotropic 
s = 1/2 Heisenberg chain (XXZ model) in the limit of high 
temperatures (/? = 0), evaluated numerically (exact diagonal- 
ization) for chain lengths L = 8, 10, . . ., 18. The curves for 
©(i) correspond to the anisotropy parameter A — 1.5 and are 
independent from the concrete choice of the field strength B. 

chain (XXZ chain) with periodic boundary conditions: 

^ = E + ^'^Ui + ^ ^l^Ui) + ^^l ■ (24) 

Here, the operators ct^ {i = x, y, z) represent the standard 
Pauli matrices (corresponding to site ^); J denotes the 
coupling strength; A is the anisotropy parameter; and B 
specifies the strength of the magnetic field. 

Fig. 1 exemplarily shows the result for A = 1.5 and 
/3 = for various lengths L ~ 8, 10, . . ., 18. Note that 
in this high temperature regime the results are indepen- 
dent from B. At short times (inset) wc find an approxi- 
mately constant I?(t) k, T) ~ 0.6 J. While the height of 
this "plateau" does not change with L, its width seems 
to increase gradually. This increase is plausible, espe- 
cially since the slope of I?(t) at large t coincides with the 
Drudc weight and is expected to vanish in the limit of 
L — > oo [2,3,6]. Remarkably, T) = 2.4 for J = 4 is very 
close to the value of 2.3 which was found in Rcf. [25] from 
a numerically involved analysis of a bath scenario. Fur- 
thermore, for a slightly different anisotropy A = 1.6 a 
similar Fig. with T> = 0.55 J is found, i.e., T> = 0.022 for 
J = 0.04. Also this result is in very good agreement with 
the value of 0.0234 in Rcf. [26] which also addresses the 
high temperature limit using bath scenarios. 

These numerical findings support the typicality state- 
ment, i.e., that the range of validity of (18) is not limited 
to the specific scenario, as explicitly analyzed in the text. 
Thus, in a forthcoming paper [27] by two of us Eq. (18) 
will be applied to various spin models in detail. 
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